Information processing apparatus, information processing method and program product for fast computation of risk measures and risk contributions

ABSTRACT

An information processing apparatus includes: a rating threshold calculating unit obtaining a probability p sk  of transition from a rating s to a rating k stored in the storage device, and calculating a rating threshold θ ir  at which the rating of a obligor i becomes a rating r or lower based on the obtained p sk , a conditional probability calculating unit obtaining a constant a i  of each obligor i stored in the storage device, and calculating a probability ps ir (z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a i , z calculated in the z calculating unit, and the θ ir  calculated in the rating threshold calculating unit, and a w ir  calculating unit calculating a difference w ir  in exposure value between adjacent ratings of the obligor i based on an exposure value V ir  of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is based upon and claims the benefit of priority of the prior Japanese Patent Application Nos. 2009-131204, filed on May 29, 2009, and 2010-095350, filed on Apr. 16, 2010, the entire contents of which are incorporated herein by reference.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to an information processing apparatus, an information processing method, and a program product.

2. Description of the Related Art

Monte Carlo methods are widely used as calculating methods of VaR (Value at Risk, a percentile according to a confidence level α) of credit risk. Further, there are DM (Default Mode) method and MTM (Mark-to-Market/market price evaluation) method as methods of modeling credit risk. The DM method is a model creating a loss distribution using only the loss in the event of default, and is used for credit risk measurement for allowing credit, including medium and small companies and retails. On the other hand, the MTM method is a model creating a profit and loss distribution considering profit and loss other than in the event of default, and is used mainly for financial products (such as bonds for example) having a marketability that requires consideration of loss due to decline in market prices.

For either model of the DM method and the MTM method, methods of calculating VaR using a Monte Carlo method are known. For example, Japanese Laid-open Patent Publication No. 2009-32237 describes explanations of modeling of the DM method and the MTM method, and a calculating method of VaR with a Monte Carlo method. Further, besides Monte Carlo methods, for example, an approach to calculate VaR analytically with respect to the DM method is disclosed in “Analytical Evaluation Method for Credit Risk of Credit Portfolio Centered on Limiting Loss Distribution and Granularity Adjustment”, IMES Discussion Paper Series 2005-J-4, Institute for Monetary and Economic Studies, Bank of Japan, July 2005.

In late years, necessity of economic capital management (operating administration in the aspect of accounting management considering viewpoints such as what degree of risk, should be taken in what sector, profitability, and the like based on risk contributions of individual companies) is discussed both domestically and globally, and importance of breaking down the total VaR into “risk contributions of individual companies relative to VaR” is increasing.

However, in the above-described Monte Carlo methods, errors occur due to random numbers when calculating risk contributions, thereby causing problems such as unstable numeric values, discrepancy between the total VaR and the sum of the contributions of individual companies, taking a long time for calculation to obtain contributions in units which are too small, requiring a high-capacity memory, and so on.

On the other hand, there are proposed techniques to calculate VaR by an analytical approach. Using the analytical approach, it is possible to calculate risk contributions of individual companies accurately and rapidly. As described above, there are the DM method and the MTM method as methods of modeling credit risk, and the IMES Discussion Paper discloses an approach to calculate VaR analytically with respect to the DM method. On the other hand, the MTM method has many parameters as compared to the DM method and provides a complicated distribution, and thus has been considered to have a difficulty to calculate VaR analytically.

Now, the above-described problems will be described in more detail using drawings. FIG. 1 is a diagram representing rating transitions and losses of a obligor i by the DM method. Reference symbol p_(i) denotes a default rate according to the rating of the obligor i. As illustrated in FIG. 1, in the DM method, ones other than default can be grouped together, and thus FIG. 1 can be expressed by a binomial distribution as illustrated in FIG. 2. FIG. 2 is a diagram illustrating an example of expressing FIG. 1 by a binomial distribution.

In the DM method, based on the rating transitions in FIG. 2, the probability distribution of loss in a portfolio is obtained as:

$L = {\sum\limits_{i}L_{i}}$ D_(i) = l_({x_(i) < θ_(i)}) L_(i) = D_(i)LGD_(i)EAD_(i) p_(i) = E[D_(i)] $l_{\{ b\}} = \left\{ \begin{matrix} 1 & {{when}\mspace{14mu} b\mspace{14mu}{is}\mspace{14mu}{true}} \\ 0 & {{when}\mspace{14mu} b\mspace{14mu}{is}\mspace{14mu}{false}} \end{matrix} \right.$ where

L is a loss of the portfolio,

L_(i) is a loss by the obligor i,

D_(i) is a random variable indicating the status of the obligor i (default or non-default),

x_(i) is a random variable indicating the enterprise value of the obligor i,

l[•] is a defined function,

p_(i) is a default rate of the obligor i,

LGD_(i) is a loss given default in exposure of the obligor i, and

EAD_(i) is an amount of exposure of the obligor i.

In practice, the default rate and the loss given default differ for every obligor included in the portfolio. Accordingly, even when an individual obligor can be expressed as in FIG. 2, the loss of the portfolio becomes a complicated probability distribution, for which it is not easy to obtain VaR.

Thus, conventionally, probability distributions and risk indicators (VaR) thereof have been calculated by simulation using a Monte Carlo method (generating random numbers corresponding to the random variable x_(i) for the number of times of trial, and calculating in the order of D_(i)→L_(i)→L).

In late years, VaR calculating methods by analytical approximation such as granularity adjustment method have been known. Accordingly, it has become possible to stably and quickly obtain not only VaR but also risk contribution of each obligor.

${{Risk}\mspace{14mu}{contribution}\mspace{14mu}{of}\mspace{14mu} a\mspace{14mu}{debtor}\mspace{14mu} i} = {{EAD}_{i}\frac{\partial{VaR}}{\partial{EAD}_{i}}}$ ${VaR} = {\sum\limits_{i}{{EAD}_{i}\frac{\partial{VaR}}{\partial{EAD}_{i}}}}$

are an expression for obtaining the risk contribution of the obligor i when obtaining VaR analytically in the DM method, and an expression for obtaining VaR analytically in the DM method.

On the other hand, in the MTM method, the diagram representing rating transitions and profits and losses of the obligor i becomes a polynomial distribution similar to FIG. 3, in which the number of parameters per obligor becomes large as compared to the DM method. Thus, there is a problem that the analytical approximation of VaR becomes quite difficult as compared to the DM method. Here, FIG. 3 is a diagram representing rating transitions and profits and losses of the obligor i by the MTM method. Reference symbol p_(ir) denotes a rating transition probability of the obligor i from the current rating to the rating r.

SUMMARY OF THE INVENTION

The present invention is made in view of such problems, and an object thereof is to allow obtaining VaR and the like quickly in analytical approximation in an MTM method.

Accordingly, an information processing apparatus of the present invention has: a z calculating unit calculating, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to α; a rating threshold calculating unit obtaining a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, and calculating a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating unit obtaining a constant a_(i) of each obligor i stored in the storage device, and calculating a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating unit, and the θ_(ir) calculated in the rating threshold calculating unit using

$\begin{matrix} {{p\;{s_{ir}(z)}} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}} & {{Expression}\mspace{14mu}(1)} \end{matrix}$ where N is a cumulative probability function of a standard normal distribution; and a w_(ir) calculating unit calculating a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w _(ir) =V _(ir) −V _(i(r+1))  Expression (2).

With the information processing apparatus having such a structure, the probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed with respect to the probability p_(sk) of transition from the rating s to the rating k can be calculated as Expression (1). Further, with the information processing apparatus having such a structure, the difference w_(ir) in exposure value between adjacent ratings of the obligor i can be calculated as Expression (2) with respect to the exposure value V_(ir) when transition of the obligor i from a current rating to the rating r or lower occurs. The point that the ps_(ir)(z) can be calculated as Expression (1) and the w_(ir) can be calculated as Expression (2) corresponds to that it is possible to express a polynomial distribution representing rating transitions and profits and losses of a obligor i in the MTM method by addition of binomial expressions. Therefore, a conditional expected value l(z) of a portfolio value V under the condition that z is fixed and a conditional variance v(z) of a portfolio value V under the condition that z is fixed can be easily (quickly) calculated, and VaR can be calculated analytically by the MTM method.

Further, an information processing method of the present invention executed by an information processing apparatus has: a z calculating step of calculating, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to α; a rating threshold calculating step of obtaining a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, and calculating a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating step of obtaining a constant a_(i) of each obligor i stored in the storage device, and calculating a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating step, and the θ_(ir) calculated in the rating threshold calculating step using

${p\;{s_{ir}(z)}} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ where N is a cumulative probability function of a standard normal distribution; and a w_(ir) calculating step of calculating a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w_(ir)=V_(ir)−V_(i(r+1)).

Further, a program product of the present invention causes a computer to function as: a z calculating unit calculating, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to α; a rating threshold calculating unit obtaining a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, and calculating a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating unit obtaining a constant a_(i) of each obligor i stored in the storage device, and calculating a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating unit, and the θ_(ir) calculated in the rating threshold calculating unit using

${p\;{s_{ir}(z)}} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ where N is a cumulative probability function of a standard normal distribution; and a w_(ir) calculating unit calculating a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w_(ir)=V_(ir)−V_(i(r+1)).

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram representing rating transitions and losses of a obligor i by a DM method;

FIG. 2 is a diagram illustrating an example of expressing FIG. 1 by a binomial distribution;

FIG. 3 is a diagram representing rating transitions and profits and losses of the obligor i by an MTM method;

FIG. 4 is a diagram illustrating an example of expressing a polynomial distribution as in FIG. 3 by addition of binomial distributions;

FIG. 5 is a diagram for explaining FIG. 4 in consideration of an enterprise value model;

FIG. 6 is a diagram illustrating an example of a hardware structure of an information processing apparatus;

FIG. 7 is a diagram illustrating an example of a software structure of the information processing apparatus;

FIG. 8 is a diagram illustrating an example of a table storing the rating transition matrix;

FIG. 9 is a diagram illustrating an example of a table storing probabilities ps_(sr) of transition from a rating s to a rating r or lower;

FIG. 10 is a diagram illustrating an example of a table storing probabilities ps_(ir) of transition of the obligor i to the rating r or lower;

FIG. 11 is a diagram illustrating an example of a table storing rating thresholds θ_(ir) at which the rating of the obligor i becomes the rating r or lower;

FIG. 12 is a diagram illustrating an example of a table storing a correlation matrix;

FIG. 13 is a diagram illustrating an example of a table storing correlation calculating information by obligor;

FIG. 14 is a diagram illustrating an example of a table storing probabilities ps_(ir)(z) that the rating of a obligor i becomes the rating r or lower under the condition that a risk factor z is fixed (conditional probability given that the risk factor z is fixed);

FIG. 15 is a diagram illustrating an example of a table storing first derivatives ps_(ir)′(z) of probabilities that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed (conditional probability given that the risk factor z is fixed);

FIG. 16 is a diagram illustrating an example of a table storing second derivatives ps_(ir)″(z) of probabilities that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed (conditional probability given that the risk factor z is fixed);

FIG. 17 is a diagram illustrating an example of a table storing current ratings of obligors i and exposure values V_(ir) when transition to the rating r occurs;

FIG. 18 is a diagram illustrating an example of a table storing differences w_(ir) in value between adjacent ratings of obligors i;

FIG. 19 is a diagram illustrating an example of a table storing conditional expected values l(z) of portfolio values V, conditional variances of the portfolio values V, and the like;

FIG. 20 is a diagram illustrating an example of a table storing expected values of exposure values by obligor, individual company risk contributions, and so on;

FIG. 21 is a flowchart illustrating an example of processing to calculate a percentile of the risk factor z corresponding to a confidence level α;

FIG. 22 is a flowchart illustrating an example of processing to calculate a rating threshold θ_(ir) at which the rating of a obligor i becomes the rating r or lower;

FIG. 23 is a flowchart illustrating an example of processing to calculate a conditional probability or the like;

FIG. 24 is a flowchart illustrating an example of processing to calculate a difference w_(ir) in value between adjacent ratings of a obligor i;

FIG. 25 is a flowchart illustrating an example of processing to calculate a percentile qα of a portfolio value or the like according to the confidence level α;

FIG. 26 is a flowchart illustrating an example of processing to calculate an expected value E[V_(i)] of an exposure value of a obligor i; and

FIG. 27 is a flowchart illustrating an example of processing to calculate an individual company risk contribution and a portfolio VaR.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Hereinafter, an embodiment of the present invention will be described based on the drawings.

Overview of this Embodiment

First, an overview of this embodiment will be described using FIG. 4 and FIG. 5. Incidentally, the processing illustrated below is executed by an information processing apparatus (computer) which will be described later.

As already described, a diagram representing rating transitions and profits and losses of a obligor i in the MTM method becomes a polynomial distribution like FIG. 3. However, when a probability ps_(ir) of transition of a obligor i from the current rating to a rating r or lower is obtained as described later corresponding to a transition probability p_(ir) of transition of the obligor i from the current rating to the rating r in a Merton model (enterprise value model), and a difference w_(ir) in exposure value between adjacent ratings of the obligor i is obtained as described later corresponding to an exposure value V_(ir) of the obligor when transition from the current rating to the rating r occurs, the polynomial distribution like FIG. 3 can be expressed by addition of binomial distributions as illustrated in FIG. 4.

Here, FIG. 5 is a diagram for explaining FIG. 4 in consideration of an enterprise value model.

In the MTM method, the information processing apparatus obtains the transition probability p_(ir) according to the rating of a obligor i from a probability transition matrix stored in a storage device such as an HD, and a rating threshold θ_(ir) is obtained so that the relation of p _(ir) =E[l _({θ) _(i(r−1)) _(≦x) _(i) _(<θ) _(ir) _(}) ]=N(θ_(ir))−N(θ_(i(r−1)) holds true, using

${p\; s_{i\; r}} = {\sum\limits_{k = 1}^{r}p_{i\; k}}$ θ_(i r) = N⁻¹(p s_(i r)), where ps_(ir) is a probability of transition of the obligor i to the rating r or lower,

N(•) is a cumulative probability function of a standard normal distribution, and

N⁻¹(•) is an inverse function of NH.

Incidentally, here the smaller the r is, the lower the rating it denotes.

Using θ_(ir), the function indicating a transition to rating BB can be represented as 42 in FIG. 5.

Here, using 43 and 44 in FIG. 5, 42 of FIG. 5 can be expressed as 43-44 of FIG. 5. That is, 43 and 44 of FIG. 5 are forms similar to D_(i) used in the DM method, and corresponds to the binomial distribution. This facilitates calculation of an expected value and a distribution. Besides the rating BB, it can be expanded as follows by expressions using S_(ir) and D_(ir).

D_(i r) = 1_({x_(i) < θ_(i r)}) p s_(i r) = E[D_(i r)] S_(i r) = D_(i r) − D_(i(r − 1)) $V_{i} = {{\sum\limits_{r}{S_{i\; r}V_{i\; r}}} = {\sum\limits_{r}{D_{i\; r}w_{i\; r}}}}$ This corresponds to expressing a polynomial distribution by addition of binomial distributions as in FIG. 4. Incidentally, mathematically a polynomial distribution cannot always be expressed by addition of binomial distributions, and the relational expressions illustrated above are relations derived by devising an expanding method in an enterprise value model.

Here, w _(ir) =V _(ir) −V _(i(r+1)),

S_(ir) is a random variable denoting the status of a obligor i (whether R_(i) is r or not),

R_(i) is a rating of the obligor i after rating transition,

D_(ir) is a random variable indicating the status of the obligor i (whether or not R_(i) is equal to or lower than rating r),

V_(i) is a value of exposure after rating transition of the obligor i, and

V_(ir) is a value of exposure when the rating of the obligor i is the rating r.

In calculation with the enterprise value model, generally an expressing method of x_(i) like the following expression for example is used. x _(i)=√{square root over (a _(i))}z+√{square root over (1−a _(i))}ζ_(i)

Here,

z is a risk factor,

a_(i) is a constant representing the dependence of an enterprise value x_(i) on z, and

ζ_(i) is a factor not expressed by z in an enterprise value x_(i) of the obligor i.

In the above expression, a conditional probability when the random variable z is fixed to a certain value as 45 in FIG. 5 is considered. Here, an expected value l(z) and a variance v(z) of the portfolio value V are needed for analytical approximation of VaR. However, expressing a polynomial distribution by addition of binomial distributions as described above is also possible here, and it is possible to obtain l(z) and v(z) by calculation based on a probability ps_(ir)(z) of transition of the obligor to the rating r or lower.

For example, in the case of multifactor with the model being expanded or in the case where it is not necessarily be Corr(ζ_(i), ζ_(j))=0, particularly calculation of v(z) and v′(z) becomes difficult, but the possibility of expression by binomial distributions facilitates calculation. Describing more specifically, it is difficult to perform the calculations as in the supplemental theories 4, 5 of the IMES Discussion Paper in the DM method as the polynomial distributions as they are similarly for the MTM method, but expression by addition of binomial distributions facilitates these calculations.

p s_(i r)(z) = E[D_(i r)❘z] ${{p\;{s_{i\; r}^{\prime}(x)}} = {\frac{\mathbb{d}}{\mathbb{d}z}p\;{s_{i\; r}(z)}}},{{p\;{s_{i\; r}^{''}(z)}} = {\frac{\mathbb{d}^{2}}{\mathbb{d}z^{2}}p\;{s_{i\; r}(z)}}}$ l(z) = E[V❘z] ${{l^{\prime}(z)} = {\frac{\mathbb{d}}{\mathbb{d}z}{l(z)}}},{{l^{''}(z)} = {\frac{\mathbb{d}^{2}}{\mathbb{d}z^{2}}{l(z)}}}$ v(z) = Var[V❘z] ${v^{\prime}(z)}\frac{\mathbb{d}}{\mathbb{d}z}{v(z)}$

As described above, the conditional expected value l(z) and the variance v(z) of the portfolio value and derivatives of them, which are required for calculating granularity adjustment, can be calculated also in the case of polynomial distributions of the MTM method. Consequently, a portfolio VaR (percentile qα according to a confidence level α) and a risk contribution of an individual company can be calculated.

$\begin{matrix} {{\Delta\;{q_{\alpha}(z)}} = {{- \frac{1}{2{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)}} \\ {{q_{\alpha} = \left( {{l(z)} + {\Delta\;{q_{\alpha}(z)}}} \right)}}_{z = {N^{- 1}{({1 - \alpha})}}} \end{matrix}$

As the analytical approximation as above, besides the above ones, saddle point approximation and the like are known. In the saddle point approximation and the like, numerical integration may be required other than the analytical method (closed-form). In this embodiment, the numerical integration is not necessary as described above, and a VaR and an individual company risk contribution can be calculated by only the closed-form. That is, according to this embodiment, it is possible to rapidly calculate the VaR and rapidly and stably calculate the individual company risk contribution.

Further, according to this embodiment, a granularity adjustment method can be applied to the MTM method. In calculation by the granularity adjustment method, dividing into a term that can be ignored and a term that cannot be ignored is possible in a limiting loss distribution (assuming that there is no credit concentration), and thus it is useful in view of grasping the influence of credit concentration.

Incidentally, when a VaR obtained by a Monte Carlo method can be obtained, this VaR may be assumed as a regular value, and the individual company risk contribution can be corrected by multiplying a correction ratio c as shown by the following expressions so that a VaR obtained by analytical approximation matches a VaR obtained by the Monte Carlo method.

c=VaR by Monte Carlo/VaR by analytical approximation Individual risk contribution of a obligor i c·(individual company risk contribution of the obligor obtained by the analytical approximation)

<Hardware Structure>

FIG. 6 is a diagram illustrating an example of a hardware structure of an information processing apparatus. As illustrated in FIG. 6, an information processing apparatus 1 includes, as hardware components, an input device 11, a display device 12, a recording medium drive device 13, a ROM (Read Only Memory) 15, a RAM (Random Access Memory) 16, a CPU (Central Processing Unit) 17, an interface device 18, and an HD (Hard Disk) 19.

The input device 11 is formed of a keyboard, a mouse, and/or the like operated by an operator (or user) of the information processing apparatus 1, and is used for inputting various operation information or the like to the information processing apparatus 1. The display device 12 is formed of a display or the like used by the user of the information processing apparatus 1, and is used for displaying various information (or screen) and the like. The interface device 18 is an interface connecting the information processing apparatus 1 to a network or the like.

A program product related to flowcharts which will be described later is provided to the information processing apparatus 1 by, for example, a recording medium 14 such as a CD-ROM, or is downloaded via a network or the like. The recording medium 14 is set to the recording medium drive device 13, and the program product is installed from the recording medium 14 to the HD 19 via the recording medium drive device 13.

The ROM 15 stores a program and the like which are read first when the power of the information processing apparatus 1 is turned on. The RAM 16 is a main memory of the information processing apparatus 1. The CPU 17 reads a program from the HD 19 as necessary, stores it in the RAM 16, and executes this program, so as to provide part of functions which will be described later, or executes a flowchart or the like which will be described later. The RAM 16 stores, for example, a table and/or the like which will be described later.

<Software Structure>

FIG. 7 is a diagram illustrating a software structure of the information processing apparatus.

As illustrated in FIG. 7, the information processing apparatus 1 includes, as software components, a z calculating unit 21, a ps_(sr) calculating unit 22, a ps_(ir) calculating unit 23, a θ_(ir) calculating unit 24, a conditional probability calculating unit 25, a w_(ir) calculating unit 26, a conditional expected value calculating unit 27, a conditional variance calculating unit 28, a qα calculating unit 29, an E[V_(i)] calculating unit 30, a risk contribution calculating unit 31, and a portfolio VaR calculating unit 32.

The z calculating unit 21 obtains a confidence level α as a parameter from the RAM 16 or the like for example, and calculates the percentile of the risk factor z corresponding to the obtained confidence level α using z=N ⁻¹(1−α). Incidentally, the parameter may be inputted by the operator of the information processing apparatus 1 via the input device 11 or the like and stored in the RAM 16 or the like, or may be stored in advance in the HD 19 or the like. The same applies to the parameters below.

The ps_(sr) calculating unit 22 obtains a rating transition matrix (probability p_(sr) of transition from the current rating s to the rating r) illustrated in FIG. 8 as a parameter from the RAM 16 or the like for example, obtains a probability ps_(sr) of transition from the current rating s to the rating r or lower based on the obtained probability p_(sr) of transition from the current rating s to the rating r, and substitutes it in the corresponding part of a table as illustrated in FIG. 9. Here, FIG. 8 is a diagram illustrating an example of a table storing the rating transition matrix. FIG. 9 is a diagram illustrating an example of a table storing probabilities ps_(sr) of transition from the rating s to the rating r or lower.

The ps_(ir) calculating unit 23 calculates the probability ps_(ir) of transition of a obligor i to the rating r or lower based on the probability ps_(sr) of transition from the rating s to the rating r or lower as illustrated in FIG. 9 calculated in the ps_(sr) calculating unit 22. Describing more specifically, the ps_(ir) calculating unit 23 obtains the current rating of the obligor i from rating information by obligor as illustrated in FIG. 17 which will be described later, obtains the corresponding probability ps_(sr) of transition from the rating s to the rating r or lower from FIG. 9 (for example, obtains the portion surrounded by a dotted line in FIG. 9) based on the obtained rating information, and substitutes it in a table as illustrated in FIG. 10 (for example, in the part surrounded by a dotted line in FIG. 9). FIG. 10 is a diagram illustrating an example of a table storing probabilities ps_(ir) of transition of the obligor i to the rating r or lower.

Here, when there are plural rating systems, there exist tables corresponding to FIG. 8 and FIG. 9 matching the number of rating systems. In this case, the ps_(ir) calculating unit 23 obtains a value with reference to FIG. 8 and FIG. 9 of the corresponding rating system for each obligor.

The ƒ_(ir) calculating unit 24 obtains the probability ps_(ir) of transition of a obligor i to the rating r or lower as illustrated in FIG. 10 calculated in the ps_(ir) calculating unit 23 from the RAM 16 or the like, calculates a rating threshold θ_(ir) at which the rating of the obligor i becomes the rating r or lower based on the obtained probability ps_(ir) using θ_(ir) =N ⁻¹(ps _(ir)) and substitutes it in the corresponding part of a table as illustrated in FIG. 11. FIG. 11 is a diagram illustrating an example of a table storing rating thresholds θ_(ir) at which the rating of the obligor i becomes the rating r or lower.

The conditional probability calculating unit 25 obtains a value a_(i) as a parameter given in units of obligors from the RAM 16 or the like for example, and calculates a probability ps_(ir)(z) that the rating of a obligor i changes to the rating r or lower under the condition that the risk factor z is fixed based on the obtained value a_(i) given in units of obligors, the percentile of the risk factor z corresponding to the confidence level α calculated in the z calculating unit 21, and θ_(ir) calculated in the θ_(ir) calculating unit 24 using

${p\;{s_{i\; r}(z)}} = {N\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 14.

Incidentally, here the value a_(i) given in units of obligors is calculated in units of obligors using a correlation matrix illustrated in FIG. 12 and correlation calculation information by obligor illustrated in FIG. 13 (a specific calculating method is shown in the IMES Discussion Paper or the like, and thus details are omitted in this embodiment).

FIG. 14 is a diagram illustrating an example of a table storing probabilities ps_(ir)(z) that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed (conditional probability given that the risk factor z is fixed).

Further, the conditional probability calculating unit 25 similarly calculates a first derivative ps_(ir)′(z) of the probability ps_(ir)(z) that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed using

${p\;{s_{i\; r}^{\prime}(z)}} = {{- \sqrt{\frac{a_{i}}{1 - a_{i}}}}{n\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 15.

FIG. 15 is a diagram illustrating an example of a table storing first derivatives ps_(ir)′(z) of probabilities that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed (conditional probability given that the risk factor z is fixed).

Further, the conditional probability calculating unit 25 similarly calculates a second derivative ps_(ir)″(z) of the probability ps_(ir)(z) that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed using

${p\;{s_{i\; r}^{''}(z)}} = {{{- \frac{a_{i}}{1 - a_{i}}} \cdot \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}}}{n\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 16.

FIG. 16 is a diagram illustrating an example of a table storing second derivatives ps_(ir)″(z) of probabilities that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed (conditional probability given that the risk factor z is fixed).

The w_(ir) calculating unit 26 calculates a difference w_(ir) in value between adjacent ratings of a obligor i based on the current rating of the obligor and the exposure value V_(ir) when transition to the rating r occurs using w _(ir) =V _(ir) −V _(i(r+1)) as illustrated in FIG. 17 as a parameter, and substitutes the calculated w_(ir) in the corresponding position of a table as illustrated in FIG. 18.

FIG. 17 is a diagram illustrating an example of a table storing current ratings of obligors i and exposure values V_(ir) when transition to the rating r occurs. FIG. 18 is a diagram illustrating an example of a table storing differences w_(ir) in value between adjacent ratings of obligors i.

The conditional expected value calculating unit 27 calculates a conditional expected value l(z) of a portfolio value V and a conditional expected value l_(i)(z) of the portfolio value V of each obligor based on the difference w_(ir) in value between adjacent ratings of a obligor i as illustrated in FIG. 18 calculated in the w_(ir) calculating unit 26, and the probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under the condition that the risk factor z is fixed as illustrated in FIG. 14 calculated in the conditional probability calculating unit 25 using

$\begin{matrix} {{{l_{i}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}(z)}}}}{{l(z)} = {{E\left\lbrack {V❘z} \right\rbrack} = {\sum\limits_{i}{l_{i}(z)}}}}} & {{Expression}\mspace{14mu} 3} \end{matrix}$ and substitutes the l_(i)(z) in the corresponding part of a table as illustrated in FIG. 19 for example.

FIG. 19 is a diagram illustrating an example of a table storing conditional expected values l(z) of portfolio values V, conditional variances of the portfolio values V, and the like.

Further, the conditional expected value calculating unit 27 similarly calculates a first derivative l′(z) of a conditional expected value l(z) of a portfolio value V and a first derivative l_(i)′(z) of a conditional expected value l_(i)(z) of the portfolio value V of each obligor using

${l_{i}^{\prime}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}^{\prime}(z)}}}$ ${l^{\prime}(z)} = {{\frac{\mathbb{d}}{\mathbb{d}z}{E\left\lbrack {V❘z} \right\rbrack}} = {\sum\limits_{i}{l_{i}^{\prime}(z)}}}$ and substitutes the l_(i)′(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

Further, the conditional expected value calculating unit 27 similarly calculates a second derivative l″(z) of a conditional expected value l(z) of a portfolio value V and a second derivative l_(i)″(z) of a conditional expected value l_(i)(z) of the portfolio value V of each obligor using

${l_{i}^{''}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}^{''}(z)}}}$ ${l^{''}(z)} = {{\frac{\mathbb{d}^{2}}{\mathbb{d}z^{2}}{E\left\lbrack {V❘z} \right\rbrack}} = {\sum\limits_{i}{l_{i}^{''}(z)}}}$ and substitutes the l_(i)″(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

The conditional variance calculating unit 28 calculates a conditional variance v(z) of a portfolio value V and a partial distribution v_(i)(z) of each obligor of v(z) based on the difference w_(ir) in value between adjacent ratings of a obligor i as illustrated in FIG. 18 calculated in the w_(ir) calculating unit 26, and the probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under the condition that the risk factor z is fixed as illustrated in FIG. 14 calculated in the conditional probability calculating unit 25 using

$\begin{matrix} {{{v_{i}(z)} = {V_{i}\frac{\partial}{\partial V_{i}}{{Var}\left\lbrack {V❘z} \right\rbrack}}}{{v(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}(z)}}}}} & {{Expression}\mspace{14mu}(4)} \end{matrix}$ and substitutes the v_(i)(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

Further, the conditional variance calculating unit 28 similarly calculates a first derivative v′(z) of a conditional variance v(z) of a portfolio value V and a partial distribution v_(i)′(z) of each obligor of v′(z) using

${v_{i}^{\prime}(z)} = {V_{i}{\frac{\partial}{\partial V_{i}} \cdot \frac{\mathbb{d}}{\mathbb{d}z}}{{Var}\left\lbrack {V❘z} \right\rbrack}}$ ${v^{\prime}(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}^{\prime}(z)}}}$ and substitutes the v_(i)′(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

The qα calculating unit 29 calculates a percentile qα or the like of a portfolio value according to the confidence level α based on the conditional expected value l_(i)(z), the first derivative l_(i)′(z), and the second derivative l_(i)″(z) of a portfolio value V of each obligor calculated in the conditional expected value calculating unit 27 and stored in the table as illustrated in FIG. 19, and the v_(i)(z) and the first derivative v_(i)′(z) calculated in the conditional variance calculating unit 28 and stored in the table as illustrated in FIG. 19 using

$\begin{matrix} {{\Delta\;{q_{\alpha\; i}(z)}} = {V_{i}\frac{\partial}{\partial V_{i}}\left\{ {{- \frac{1}{2{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)} \right\}}} \\ {{\Delta\;{q_{\alpha}(z)}} = {\sum\limits_{i}{\Delta\; q_{\alpha\; i}}}} \\ {{q_{\alpha\; i} = \left( {{l_{i}(z)} + {\Delta\;{q_{\alpha\; i}(z)}}} \right)}}_{z = {N^{- 1}{({1 - \alpha})}}} \\ {q_{\alpha} = {\sum\limits_{i}q_{\alpha\; i}}} \end{matrix}$ and substitutes the calculated qα_(i) or Δqα_(i)(z) in the corresponding part of the table as illustrated in FIG. 19.

The E[V_(i)] calculating unit 30 obtains a rating transition matrix (probability p_(sr) of transition from the current rating s to the rating r) illustrated in FIG. 8 as a parameter from the RAM 16 or the like for example and obtains an exposure value V_(ir) of a obligor when transition from the current rating to the rating r occurs as illustrated in FIG. 17 as a parameter from the RAM 16 or the like for example, and calculates an expected value E[V_(i)] of an exposure value of a obligor i based on the obtained probability p_(sr) of transition from the current rating to the rating r and the obtained current rating of the obligor i and exposure value V_(ir) when transition to the rating r occurs using

${E\left\lbrack V_{i} \right\rbrack} = {\sum\limits_{r}{p_{i\; r}V_{i\; r}}}$ and substitutes the calculated E[V_(i)] in the corresponding position of a table as illustrated in FIG. 20.

FIG. 20 is a diagram illustrating an example of a table storing expected values of exposure values by obligor, individual company risk contributions, and so on.

The risk contribution calculating unit 31 calculates a risk contribution (individual company risk contribution) of a obligor i based on the E[V_(i)] calculated in the E[V_(i)] calculating unit 30 and stored in the table as illustrated in FIG. 20 and the qα_(i) calculated in the qα calculating unit 29 and stored in the table as illustrated in FIG. 19 using risk contribution of a obligor i=qα _(i) −E[V _(i)] for example and substitutes the calculated risk contribution of the obligor i in the corresponding part of the table as illustrated in FIG. 20.

The portfolio VaR calculating unit 32 calculates a portfolio VaR based on the risk contribution of a obligor i calculated in the risk contribution calculating unit 31 and stored in the table as illustrated in FIG. 20 using

${{Portfolio}\mspace{14mu}{VaR}} = {{{q\;\alpha} - {\sum\limits_{i}{E\left\lbrack V_{i} \right\rbrack}}} = {\sum\limits_{i}\mspace{14mu}{{risk}\mspace{14mu}{contribution}\mspace{14mu}{of}\mspace{14mu}{an}\mspace{14mu}{obligor}\mspace{14mu} i}}}$ for example and substitutes the calculated risk contribution of the obligor i in the corresponding part of the table as illustrated in FIG. 20.

Hereinafter, processing of the information processing apparatus will be described using flowcharts.

FIG. 21 is a flowchart illustrating an example of processing to calculate the percentile of the risk factor z corresponding to the confidence level α.

In step S10, the z calculating unit 21 obtains a confidence level α as a parameter from the RAM 16 or the like for example, and calculates the percentile of the risk factor z corresponding to the obtained confidence level α using z=N ⁻¹(1−α).

FIG. 22 is a flowchart illustrating an example of processing to calculate the rating threshold θ_(ir) at which the rating of a obligor i becomes the rating r or lower.

In step S20, the ps_(sr) calculating unit 22 obtains a rating transition matrix (probability p_(sr) of transition from the current rating s to the rating r) illustrated in FIG. 8 as a parameter from the RAM 16 or the like for example, obtains a probability ps_(sr) of transition from the current rating s to the rating r or lower based on the obtained probability p_(sr) of transition from the current rating s to the rating r, and substitutes it in the corresponding part of a table as illustrated in FIG. 9.

In step S21, the ps_(ir) calculating unit 23 calculates the probability ps_(ir) of transition of a obligor i to the rating r or lower based on the probability ps_(sr) of transition from the rating s to the rating r or lower as illustrated in FIG. 9 calculated in step S20, and substitutes it in the corresponding part of a table as illustrated in FIG. 10.

In step S22, the θ_(ir) calculating unit 24 obtains the probability ps_(ir) of transition of a obligor i to the rating r or lower as illustrated in FIG. 10 calculated in step S21 from the RAM 16 or the like, calculates a rating threshold θ_(ir) at which the rating of the obligor i becomes the rating r or lower based on the obtained probability ps_(ir) using θ_(ir) =N ⁻¹(ps _(ir)) and substitutes it in the corresponding part of a table as illustrated in FIG. 11.

FIG. 23 is a flowchart illustrating an example of processing to calculate the conditional probability or the like.

In step S30, the conditional probability calculating unit 25 obtains a value a_(i) as a parameter given in units of obligors from the RAM 16 or the like for example, and calculates a probability ps_(ir)(z) that the rating of a obligor i changes to the rating r or lower under the condition that the risk factor z is fixed based on the obtained a_(i), the percentile of the risk factor z corresponding to the confidence level α calculated in the z calculating unit 21, and θ_(ir) calculated in the θ_(ir) calculating unit 24 using

${p\;{s_{i\; r}(z)}} = {N\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 14.

In step S31, the conditional probability calculating unit 25 similarly calculates a first derivative ps_(ir)′(z) of the probability ps_(ir)(z) that the rating of a obligor i becomes the rating r or lower under the condition that the risk factor z is fixed using

${p\;{s_{i\; r}^{\prime}(z)}} = {{- \sqrt{\frac{a_{i}}{1 - a_{i}}}}{n\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 15.

In step S32, the conditional probability calculating unit 25 similarly calculates a second derivative ps_(ir)″(z) of the probability ps_(ir)(z) that the rating of a obligor I becomes the rating r or lower under the condition that the risk factor z is fixed using

${p\;{s_{i\; r}^{''}(z)}} = {{{- \frac{a_{i}}{1 - a_{i}}} \cdot \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}}}{n\left( \frac{\theta_{i\; r} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}}$ and substitutes it in the corresponding part of a table as illustrated in FIG. 16.

FIG. 24 is a flowchart illustrating an example of processing to calculate the difference w_(ir) in value between adjacent ratings of a obligor i.

In step S40, the w_(ir) calculating unit 26 calculates a difference w_(ir) in value between adjacent ratings of a obligor I based on the current rating of the obligor I and the exposure value V_(ir) when transition to the rating r occurs using w _(ir) =V _(ir) −V _(i(r+1)) as illustrated in FIG. 17 as a parameter, and substitutes the calculated w_(ir) in the corresponding position of a table as illustrated in FIG. 18.

FIG. 25 is a flowchart illustrating an example of processing to calculate the percentile qα of a portfolio value or the like according to the confidence level α.

In step S50, the conditional expected value calculating unit 27 calculates a conditional expected value l(z) of a portfolio value V and a conditional expected value l_(i)(z) of the portfolio value V of each obligor based on the difference w_(ir) in value between adjacent ratings of a obligor I as illustrated in FIG. 18 calculated in the processing of FIG. 24, and the probability ps_(ir)(z) that the rating of the obligor I becomes the rating r or lower under the condition that the risk factor z is fixed as illustrated in FIG. 14 calculated in the processing of FIG. 23 using

${l_{i}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}(z)}}}$ ${l(z)} = {{E\left\lbrack {V❘z} \right\rbrack} = {\sum\limits_{i}{l_{i}(z)}}}$ and substitutes the l_(i)(z) in the corresponding part of a table as illustrated in FIG. 19 for example.

Further, the conditional expected value calculating unit 27 similarly calculates a first derivative l′(z) of a conditional expected value l(z) of a portfolio value V and a first derivative l_(i)′(z) of a conditional expected value l_(i)(z) of the portfolio value V of each obligor using

${l_{i}^{\prime}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}^{\prime}(z)}}}$ ${l^{\prime}(z)} = {{\frac{\mathbb{d}}{\mathbb{d}z}{E\left\lbrack {V❘z} \right\rbrack}} = {\sum\limits_{i}{l_{i}^{\prime}(z)}}}$ and substitutes the l_(i)′(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

Further, the conditional expected value calculating unit 27 similarly calculates a second derivative l″(z) of a conditional expected value l(z) of a portfolio value V and a second derivative l_(i)″(z) of a conditional expected value l_(i)(z) of the portfolio value V of each obligor using

${l_{i}^{''}(z)} = {\sum\limits_{r}{w_{i\; r}p\;{s_{i\; r}^{''}(z)}}}$ ${l^{''}(z)} = {{\frac{\mathbb{d}^{2}}{\mathbb{d}z^{2}}{E\left\lbrack {V❘z} \right\rbrack}} = {\sum\limits_{i}{l_{i}^{''}(z)}}}$ and substitutes the l_(i)″(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

In step S51, the conditional variance calculating unit 28 calculates a conditional variance v(z) of a portfolio value V and a partial distribution v_(i)(z) of each obligor of v(z) based on the difference w_(ir) in value between adjacent ratings of a obligor i as illustrated in FIG. 18 calculated in the processing of FIG. 24, and the probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under the condition that the risk factor z is fixed as illustrated in FIG. 14 calculated in the processing of FIG. 23 using

${v_{i}(z)} = {V_{i}\frac{\partial}{\partial V_{i}}{{Var}\left\lbrack {V❘z} \right\rbrack}}$ ${v(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}(z)}}}$ and substitutes the v_(i)(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

Further, the conditional variance calculating unit 28 similarly calculates a first derivative v′(z) of a conditional variance v(z) of a portfolio value V and a partial distribution v_(i)′(z) of each obligor of v′(z) using

${v_{i}^{\prime}(z)} = {V_{i}{\frac{\partial}{\partial V_{i}} \cdot \frac{\mathbb{d}}{\mathbb{d}z}}{{Var}\left\lbrack {V❘z} \right\rbrack}}$ ${v^{\prime}(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}^{\prime}(z)}}}$ and substitutes the v_(i)′(z) in the corresponding part of the table as illustrated in FIG. 19 for example.

Here, the order of the processing of step S50 and the processing of step S51 is not in question. Step S51 may be earlier than step S50, or may be later than step S50. Further, they may be performed simultaneously.

In step S52, the qα calculating unit 29 calculates a percentile qα or the like of a portfolio value according to the confidence level α based on the conditional expected value l_(i)(z), the first derivative l_(i)′(z), and the second derivative l_(i)″(z) of a portfolio value V of each obligor calculated in step S50 and stored in the table as illustrated in FIG. 19, and the v_(i)(z) and the first derivative v_(i)′(z) calculated in step S51 and stored in the table as illustrated in FIG. 19 using

$\begin{matrix} {{\Delta\;{q_{\alpha\; i}(z)}} = {V_{i}\frac{\partial}{\partial V_{i}}\left\{ {{- \frac{1}{2{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)} \right\}}} \\ {{\Delta\;{q_{\alpha}(z)}} = {\sum\limits_{i}{\Delta\; q_{\alpha\; i}}}} \\ {{q_{\alpha\; i} = \left( {{l_{i}(z)} + {\Delta\;{q_{\alpha\; i}(z)}}} \right)}}_{z = {N^{- 1}{({1 - \alpha})}}} \\ {q_{\alpha} = {\sum\limits_{i}q_{\alpha\; i}}} \end{matrix}$ and substitutes the calculated qα_(i) or Δaα_(i)(z) in the corresponding part of the table as illustrated in FIG. 19.

FIG. 26 is a flowchart illustrating an example of processing to calculate the expected value E[V_(i)] of the exposure value of a obligor i.

In step S60, the E[V_(i)] calculating unit 30 obtains a rating transition matrix (probability p_(sr) of transition from the current rating s to the rating r) illustrated in FIG. 8 as a parameter from the RAM 16 or the like for example and obtains an exposure value V_(ir) of a obligor i when transition from the current rating to the rating r occurs as illustrated in FIG. 17 as a parameter from the RAM 16 or the like for example, and calculates an expected value E[V_(i)] of an exposure value of a obligor i based on the obtained probability p_(sr) of transition from the current rating to the rating r and the obtained current rating of the obligor i and exposure value V_(ir) when transition to the rating r occurs using

${E\left\lbrack V_{i} \right\rbrack} = {\sum\limits_{r}{p_{i\; r}V_{i\; r}}}$ and substitutes the calculated E[V_(i)] in the corresponding position of a table as illustrated in FIG. 20.

FIG. 27 is a flowchart illustrating an example of processing to calculate the individual company risk contribution and the portfolio VaR.

In step S70, the risk contribution calculating unit 31 calculates a risk contribution (individual company risk contribution) of a obligor i based on the E[V_(i)] calculated in the processing of FIG. 26 and stored in the table as illustrated in FIG. 20 and the qα_(i) calculated in the processing of FIG. 25 and stored in the table as illustrated in FIG. 19 using risk contribution of a obligor i=qα _(i) −E[V _(i)] for example and substitutes the calculated risk contribution of the obligor i in the corresponding part of the table as illustrated in FIG. 20.

In step S71, the portfolio VaR calculating unit 32 calculates a portfolio VaR based on the risk contribution of a obligor i calculated in step S70 and stored in the table as illustrated in FIG. 20 using

${{Portfolio}\mspace{14mu}{VaR}} = {{{q\;\alpha} - {\sum\limits_{i}{E\left\lbrack V_{i} \right\rbrack}}} = {\sum\limits_{i}\mspace{14mu}{{risk}\mspace{14mu}{contribution}\mspace{14mu}{of}\mspace{14mu}{an}\mspace{14mu}{obligor}\mspace{14mu} i}}}$ for example and substitutes the calculated risk contribution of the obligor i in the corresponding part of the table as illustrated in FIG. 20.

As described above, the information processing apparatus of this embodiment has a technical characteristic in that it calculates analytical approximation of VaR of the portfolio value V in the MTM method using the difference w_(ir) in exposure value between adjacent ratings. Hereinafter, this point will be described in more detail.

First, calculation of VaR by the Monte Carlo method in each of the DM method and the MTM method will be described.

Expressions for performing a simulation by the Monte Carlo method in the DM method are shown below.

$\begin{matrix} {{\theta_{i} = {N^{- 1}\left( p_{i} \right)}}{D_{i} = 1_{\{{x_{i} < \theta_{i}}\}}}{1_{\{ b\}} = \left\{ {{\begin{matrix} 1 & {{when}\mspace{14mu} b\mspace{14mu}{is}\mspace{14mu}{true}} \\ 0 & {{when}\mspace{14mu} b\mspace{14mu}{is}\mspace{14mu}{false}} \end{matrix}L_{i}} = {{D_{i}{LGD}_{i}{EAD}_{i}L} = {\sum\limits_{i}L_{i}}}} \right.}} & {{Expression}\mspace{14mu}(5)} \end{matrix}$

Here,

L is a loss of the portfolio,

L_(i) is a loss by the obligor i,

D_(i) is a random variable indicating the status of the obligor i (default or non-default),

x_(i) is a random variable indicating the enterprise value of the obligor i,

l[•] is a defined function,

x_(i)<θ_(i) is whether the enterprise value x_(i) of a obligor i becomes lower than the default threshold θ_(i) or not

LGD_(i) is a loss given default in exposure of the obligor i, and

EAD_(i) is an amount of exposure of the obligor i,

p_(i) is a default rate of the obligor i

N(•) is a cumulative probability function of a standard normal distribution, and

N⁻¹(•) is an inverse function of N(•).

Expressions for performing a simulation by the Monte Carlo method in the MTM method are shown below.

$\begin{matrix} {{{p\; s_{i\; r}} = {\sum\limits_{k = 1}^{r}p_{i\; k}}}{\theta_{i\; r} = {N^{- 1}\left( {p\; s_{i\; r}} \right)}}{S_{i\; r} = {1_{\{{R_{i} = r}\}} = {{1_{\{{\theta_{i{({r - 1})}} \leq x_{i} < \theta_{i\; r}}\}}V_{i}} = {{\sum\limits_{r}{S_{i\; r}V_{i\; r}V}} = {\sum\limits_{i}V_{i}}}}}}} & {{Expression}\mspace{14mu}(6)} \end{matrix}$

Here,

V is the value of a portfolio,

V_(i) is a value of exposure after rating transition of a obligor i

V_(ir) is a value of exposure when the obligor i has a rating r,

R_(i) is the rating of the obligor i after rating transition,

S_(ir) is a random variable indicating the status of the obligor i (whether Ri is r or not),

x_(i) is a random variable indicating the enterprise value of the obligor i,

p_(ik) is a probability of transition of the obligor to a rating k,

ps_(ir) is a probability of transition of the obligor to a rating r,

N(•) is a cumulative probability function of a standard normal distribution, and

N⁻¹(•) is an inverse function of N(•).

Note that the smaller the r is, the lower the rating it denotes in the above-described expressions.

As shown in Expression (5), in the DM method, the information processing apparatus generates random numbers corresponding to the random variable x_(i) for the number of times of trials and calculates in the order of D_(i)→L_(i)→L, to thereby calculate a probability distribution and a risk indicator (namely, VaR) thereof. Further, as shown in Expression (6), in the MTM method, the information processing apparatus generates random numbers corresponding to the random variable x_(i) for the number of times of trials and calculates in the order of S_(ir)→V_(i)→V, to thereby calculate a probability distribution and a risk indicator (namely, VaR) thereof.

Next, analytical approximation of VaR will be described. For example, the IMES Discussion Paper discloses an approach of analytical approximation of VaR in the DM method. The IMES Discussion Paper describes approximation of a quantile qα(L) of a loss ratio distribution L with qα(E[L|Z])+Δqα(L). Incidentally, the IMES Discussion Paper describes as theory 2 that an α quantile qα(E[L|Z]) of a limiting loss distribution E[L|Z] is calculated by a conditional expected loss given a 1−a quantile of Z. Further, the IMES Discussion Paper describes as theory 3 a granularity adjustment term which is a difference Δqα(L) (=qα(L)−qα(E[L|Z])) between a true loss quantile and an a quantile of a limiting loss distribution. Then the IMES Discussion Paper describes that this granularity adjustment term Δqα(L) is represented as the following expression. Here, l(z)=E[L|Z=z], v(z)=VaR[L|Z=z].

${{{\Delta\;{q_{\alpha}(L)}} = {{- \frac{1}{2{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)}}}_{z = {N^{- 1}{({1 - \alpha})}}}$

Then the IMES Discussion Paper describes that a VaR is calculated as the quantile qα(L) of the loss ratio distribution L by calculating the above expression.

Incidentally, in the Monte Carlo method, the information processing apparatus uses

$V_{i} = {\sum\limits_{r}{S_{ir}V_{ir}}}$ in the MTM method corresponding to the calculation of L_(i)=D_(i)LGD_(i)EAD_(i) of the DM method as shown in Expression (6). Accordingly, the information processing apparatus calculates the above-described granularity adjustment term using this expression in the analytical approximation of the MTM method.

Here, l(z)=E[V|z], v(z)=VaR[V|z]. This calculation of VaR[V|z] is represented more specifically as the following expression.

$\begin{matrix} {{{Var}\left\lbrack V \middle| z \right\rbrack} = {\sum\limits_{i}{\sum\limits_{j}{{Cov}\left\lbrack {V_{i},\left. V_{j} \middle| z \right.} \right\rbrack}}}} \\ {= {\sum\limits_{i}{\sum\limits_{j}{\sum\limits_{r}{\sum\limits_{s}{{Corr}\left\lbrack {{S_{ir}V_{ir}},\left. {S_{js}V_{js}} \middle| z \right.} \right\rbrack}}}}}} \\ {\sqrt{{{Var}\left\lbrack {S_{ir}V_{ir}} \middle| z \right\rbrack} \cdot {{Var}\left\lbrack {S_{js}V_{js}} \middle| z \right\rbrack}}} \end{matrix}$

Here, the amount of calculation of the part Corr, [S_(ir)V_(ir), S_(js)V_(j)|z] is proportional to the “square of (the number of obligors×the number of ratings)”. For example, when the number of obligors=10000 and the number of ratings=10, the calculation of correlation is needed to be performed the number of times equal to the square of (10000×10), that is, the tenth power of 10.

On the other hand, the information processing apparatus according to this embodiment calculates the a quantile qα of a portfolio value V using the difference w_(ir) in exposure value between adjacent ratings. At this time, the information processing apparatus calculates E[V|z] in Expression (3) or the like and VaR[V|z] in Expression (4) or the like for calculating Δqα, but also in this calculation, the amount of calculation in proportion to “the square of (the number of obligors×the number of ratings)” is needed. However, as described using FIG. 15, FIG. 16, and the like of this embodiment, the information processing apparatus calculates N⁻¹(ps_(ir)(z)) and the like besides w_(ir), ps_(ir)(z), and ps′_(ir)(z) for “the number of obligors×the number of ratings”, and stores the calculated values of correlations ρij for the amount of “the square of the number of obligors” as parameters in a storage device such as a memory or a cache. In this point, the information processing apparatus can read and reuse data of parameters from the storage device such as a memory when executing calculation of the amount in proportion to “the square of (the number of obligors×the number of ratings)” requiring these parameters. Accordingly, calculation efficiency increases significantly. Further, the information processing apparatus keeps frequently used data stored in a cache having much higher reading and writing speed than that of the memory. Thus, when these parameters are stored in the cache and the operation is increased in speed particularly using a function called SIMD (Single Instruction Multiple Data) of the CPU or the like, it is possible to prevent the reading speed of the memory from becoming a bottle neck, and reduce the idling time of the CPU or the like. That is, when the analytical approximation of VaR in the MTM method is obtained, the difference w_(ir) in exposure value between adjacent ratings can be used to increase calculation efficiency, and thus the VaR in the MTM method can be calculated rapidly.

Incidentally, v_(i)(z) and vi′(z) are represented by the following expressions.

$\begin{matrix} {{{v_{i}(z)} = {{v_{\;{\infty\; i}}(z)} + {v_{G\; i}\;(z)}}}{{v_{i}^{\prime}(z)} = {{v_{\infty\; i}^{\prime}(z)} + {v_{Gi}^{\prime}\;(z)}}}} & {{Expresssion}\mspace{14mu}(7)} \end{matrix}$

Here, v_(∞i) and v′_(∞i) are represented by the following expressions.

$\begin{matrix} {{{v_{\infty\; i}(z)} = {2{\sum\limits_{j}{\sum\limits_{r}{\sum\limits_{s}{w_{ir}{w_{js}\left( {{N_{2}\left( {{N^{- 1}\left( {p\;{s_{ir}(z)}} \right)},{N^{- 1}\left( {p\;{s_{js}(z)}} \right)},\rho_{ij}} \right)} - {p\;{s_{ir}(z)}p\;{s_{js}(z)}}} \right)}}}}}}}{{v_{\infty\; i}^{\prime}(z)} = {{2{\sum\limits_{j}{\sum\limits_{r}{\sum\limits_{s}{w_{ir}w_{js}p\;{s_{ir}^{\prime}(z)}\left( {{N\left( \frac{{N^{- 1}\left( {p\;{s_{js}(z)}} \right)} - {\rho_{ij}{N^{- 1}\left( {p\;{s_{ir}(z)}} \right)}}}{\sqrt{1 - \left( \rho_{ij} \right)^{2}}} \right)} - {p\;{s_{js}(z)}}} \right)}}}}} + {2{\sum\limits_{j}{\sum\limits_{r}{\sum\limits_{s}{w_{ir}w_{js}p\;{s_{js}^{\prime}(z)}\left( {{N\left( \frac{{N^{- 1}\left( {p\;{s_{ir}(z)}} \right)} - {\rho_{ij}{N^{- 1}\left( {p\;{s_{js}(z)}} \right)}}}{\sqrt{1 - \left( \rho_{ij} \right)^{2}}} \right)} - {p\;{s_{ir}(z)}}} \right)}}}}}}}} & {{Expression}\mspace{14mu}(8)} \end{matrix}$

Further, v_(Gi) and v′_(Gi) are represented by the following expressions.

$\begin{matrix} {{v_{Gi}(z)} = {2{\sum\limits_{r}{\sum\limits_{s}{w_{ir}{w_{is}\left( {{p\;{s_{it}(z)}} - {N_{2}\left( {{N^{- 1}\left( {p\;{s_{ir}(z)}} \right)},{N^{- 1}\left( {p\;{s_{is}(z)}} \right)},\rho_{ii}} \right)}} \right)}}}}}} & {{Expression}\mspace{14mu}(9)} \\ {{v_{Gi}^{\prime}(z)} = {2{\sum\limits_{r}{\sum\limits_{s}{w_{ir}{w_{is}\left( {{p\;{s_{it}^{\prime}(z)}} - {2p\;{s_{ir}^{\prime}(z)}{N\left( \frac{{N^{- 1}\left( {p\;{s_{is}(z)}} \right)} - {\rho_{ii}{N^{- 1}\left( {p\;{s_{ir}(z)}} \right)}}}{\sqrt{1 - \left( \rho_{ii} \right)^{2}}} \right)}}} \right)}}}}}} & \; \end{matrix}$

Parts of N₂(x, y, ρ) in Expression (8) and Expression (9) are cumulative probability functions of a bivariate standard normal distribution. Further, t=min(r, s) (one with a lower rating out of r and s). Further, ρ_(ij)=corr(x_(i), x_(j)|z) (the part ε_(i) corresponding to the idiosyncratic risk of a obligor in x_(i) is calculated as corr(ε_(i), ε_(j))=0).

Moreover, Δaα_(i)(z) is represented by, more specifically, the following expression.

$\begin{matrix} {{\Delta\;{q_{\alpha\; i}(z)}} = {{V_{i}\frac{\partial}{{\partial V_{i}}\;}\Delta\;{q_{\alpha}(z)}} = {{\frac{l_{i}^{\prime}(z)}{2\left( {l^{\prime}(z)} \right)^{2}}{v^{\prime}(z)}} - {\frac{1}{2{l^{\prime}(z)}}{v_{i}^{\prime}(z)}} - {\frac{l_{i}^{\prime}(z)}{2\left( {l^{\prime}(z)} \right)^{2}}{v(z)}z} + {\frac{1}{2{l^{\prime}(z)}}{v_{i}(z)}z} - {\frac{l_{i\;}^{\prime}(z)}{\left( {l^{\prime}(z)} \right)^{3}}{v(z)}{l^{''}(z)}} + {\frac{1}{2\left( {l^{\prime}(z)} \right)^{2}}{v_{i}(z)}{l^{''}(z)}} + {\frac{1}{2\left( {l^{\prime}(z)} \right)^{2}}{v(z)}{l_{i}^{''}(z)}}}}} & {{Expression}\mspace{14mu}(10)} \end{matrix}$

With the above-described structure, the information processing apparatus according to this embodiment uses the difference w_(ir) in exposure value between adjacent ratings so as to store N⁻¹(ps_(ir)(z)), ρij, and the like besides w_(ir), ps_(ir)(z), and ps′_(ir)(z) as parameters in the cache, which can be reused repeatedly when performing the calculation of Expression (7) to Expression (10). Thus the calculation efficiency can be increased. Accordingly, the analytical approximation of VaR in the MTM method can be calculated rapidly.

As has been described above, according to this embodiment, it becomes possible to quickly calculate the conditional expected value l(z) and the variance v(z) of a portfolio value needed for calculating granularity adjustment, and the derivatives of them for polynomial distributions in the MTM method. Further consequently, the portfolio VaR (percentile qα according to the confidence level α) and risk contributions of individual companies can be calculated rapidly.

In short, according to this embodiment, a polynomial distribution representing rating transitions and profits and losses of a obligor i in the MTM method can be expressed by addition of binomial distributions, and thereby VaR and the like can be obtained quickly by analytical approximation even in the MTM method.

Note that the processing of this embodiment can be applied to calculation of a risk indicator similar to VaR, such as an expected shortfall, or the like.

Further, the above-described l_(i)(z), l_(i)′(z), l_(i)″(z), v_(i)(z), v_(i)′(z), and qα_(i) are repeating of the same calculation pattern for each obligor. Accordingly, the conditional expected value calculating unit 27, the conditional variance calculating unit 28, the qα calculating unit 29, and the like can use the function called SIMD (Single Instruction Multiple Data) mounted in the CPU 17 or the like (hereinafter referred to as a SIMD function) to calculate the l_(i)(z), l_(i)′(z), l_(i)″(z), v_(i)(z), v_(i)′(z), and qα_(i) rapidly.

In the foregoing, the preferred embodiment of the present invention has been described in detail, but the present invention is not limited to such a specific embodiment, and can be modified and changed in various ways within the scope of the spirit of the invention described in the claims that follow.

According to the present invention, it is possible to allow obtaining VaR and the like quickly in analytical approximation in an MTM method.

The present embodiments are to be considered in all respects as illustrative and no restrictive, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. The invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. 

The invention claimed is:
 1. An information processing apparatus, comprising: a computer processor; a z calculating unit configured to calculate, via the computer processor, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to α using z=N ⁻¹(1−α); a rating threshold calculating unit configured to obtain, via the processor, a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, wherein the rating threshold calculating unit is further configured to calculate a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating unit configured to obtain, via the computer processor, a constant a_(i) of each obligor i stored in the storage device, wherein the conditional probability calculating unit is further configured to calculate a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating unit, and the θ_(ir) calculated in the rating threshold calculating unit using ${p\;{s_{ir}(z)}} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ where N is a cumulative probability function of a standard normal distribution; a w_(ir) calculating unit configured to calculate, via the computer processor, a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w_(ir)=V_(ir)−V_(i(r+1)); a conditional expected value calculating unit configured to calculate, via the computer processor, a conditional expected value 1(z) of a portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating unit and ps_(ir)(z) calculated in the conditional probability calculating unit using ${l_{i}(z)} = {\sum\limits_{r}{w_{ir}p\;{s_{ir}(z)}}}$ ${{l(z)} = {{E\left\lbrack V \middle| z \right\rbrack} = {\sum\limits_{i}{l_{i}(z)}}}};$ a conditional variance calculating unit configured to calculate, via the computer processor, a conditional variance v(z) of the portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating unit and ps_(ir)(z) calculated in the conditional probability calculating unit using ${v_{i}(z)} = {V_{i}\frac{\partial\;}{\partial V_{i}}{{Var}\left\lbrack V \middle| z \right\rbrack}}$ ${{v(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}(z)}}}};$ a qα calculating unit configured to calculate, via the computer processor, a percentile qα of a portfolio value based on a first derivative 1′(z) of 1(z) calculated in the conditional expected value calculating unit, a second derivative 1″(z) of 1(z) calculated in the conditional expected value calculating unit, v(z) calculated in the conditional variance calculating unit, and a first derivative v′(z) of v(z) calculated in the conditional variance calculating unit using ${\Delta\;{q_{\alpha\; i}(z)}} = {V_{i}\frac{\partial\;}{\partial V_{i}}\left\{ {{- \frac{1}{2{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)} \right\}}$ ${\Delta\;{q_{\alpha}(z)}} = {\sum\limits_{i}{\Delta\; q_{\alpha\; i}}}$ ${q_{\alpha\; i} = {\left. \left( {{l_{i}(z)} + {\Delta\;{q_{\alpha\; i}(z)}}} \right) \middle| {}_{z = {N^{- 1}{({1 - \alpha})}}}q_{\alpha} \right. = {\sum\limits_{i}q_{\alpha\; i}}}};$ an E[V_(i)] calculating unit configured to obtain, via the computer processor, an exposure value Vir of the obligor i when transition from a current rating to the rating r occurs, wherein the E[V_(i)] calculating unit is further configured to obtain a probability p_(ir) of transition from a rating i to the rating r stored in the storage device, and wherein the E[V_(i)] calculating unit is further configured to calculate an expected value E[V_(i)] of an exposure value of the obligor i based on the obtained V_(ir) and probability p_(ir) using ${{E\left\lbrack V_{i\;} \right\rbrack} = {\sum\limits_{r}{p_{ir}V_{ir}}}};\mspace{14mu}{and}$ a portfolio VaR calculating unit configured to calculate, via the computer processor, a portfolio VaR from the percentile q_(αi) and the expected value E[V_(i)] of the exposure value of the obligor I, wherein the rating threshold calculating unit includes: a ps_(sr) calculating unit configured to calculate, via the computer processor, a probability ps_(sr) of transition from the rating s to the rating r or lower based on the probability p_(sk), using ${ps}_{sr} = {\sum\limits_{k = 1}^{r}p_{sk}}$ a ps_(ir) calculating unit configured to obtain, via the computer processor, a probability psi, of transition of the obligor i to the rating r or lower based on the probability ps_(sr) calculated in the ps_(sr) calculating unit; and a θ_(ir) calculating unit configured to calculate, via the computer processor, the rating threshold θ_(ir) at which the rating of the obligor i becomes the rating r or lower based on the probability ps_(ir) obtained in the ps_(ir) calculating unit using θ_(ir)=N⁻¹(ps_(ir)).
 2. The information processing apparatus according to claim 1, wherein the portfolio VaR calculating unit is configured to calculate, via the computer processor, a portfolio VaR based on E[V_(i)] calculated in the E[V_(i)] calculating unit and qα calculated in the qα calculating unit using ${{Portfolio}\mspace{14mu}{VaR}} = {{q\;\alpha}\; - {\sum\limits_{i}{{E\left\lbrack V_{i} \right\rbrack}.}}}$
 3. The information processing apparatus according to claim 1, further comprising: a risk contribution calculating unit configured to calculate, via the computer processor, a risk contribution of a obligor i based on E[V_(i)] calculated in the E[V_(i)] calculating unit and qα_(i) calculated when calculating qα in the qα calculating unit using risk contribution of a obligor i=qα _(i) −E[V _(i)].
 4. The information processing apparatus according to claim 3, further comprising: a portfolio VaR calculating unit configured to calculate, via the computer processor, a portfolio VaR based on the risk contribution of the obligor i calculated in the risk contribution calculating unit using Portfolio VaR=Σrisk contribution of an obligor i.
 5. The information processing apparatus according to claim 1, wherein the conditional expected value calculating unit is configured to calculate the conditional expected value 1(z) of the portfolio value V using a SIMD function provided by a CPU, wherein the conditional variance calculating unit is configured to calculate the conditional variance v(z) of the portfolio V using the SIMD function provided by the CPU, and wherein the qα calculating unit is configured to calculate qα_(i) using the SIMD function provided by the CPU.
 6. An information processing method executed by an information processing apparatus, the method comprising: a z calculating step of calculating via a computing device, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to a using z=N ⁻¹(1−α); a rating threshold calculating step of obtaining via the computing device a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, and calculating via the computing device a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating step of obtaining a constant a_(i) of each obligor stored in the storage device, and calculating via the computing device a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating step, and the θ_(ir) calculated in the rating threshold calculating step using ${{ps}_{ir}(z)} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ where N is a cumulative probability function of a standard normal distribution; a w_(ir) calculating step of calculating via the computing device a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w_(ir)=V_(ir)−V_(i(r+1)); a conditional expected value calculating step of calculating, via the computing device, a conditional expected value 1(z) of a portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating step and ps_(ir)(z) calculated in the conditional probability calculating step using $\begin{matrix} {{l_{i}(z)} = {\sum\limits_{r}{w_{ir}p\;{s_{ir}(z)}}}} \\ {{{l(z)} = {{E\left\lbrack V \middle| z \right\rbrack} = {\sum\limits_{i}{l_{i}(z)}}}};} \end{matrix}$ a conditional variance calculating step of calculating, via the computing device, a conditional variance v(z) of the portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating step and ps_(ir)(z) calculated in the conditional probability calculating step using $\begin{matrix} {{v_{i}(z)} = {V_{i}\frac{\partial\;}{\partial V_{i}}{{Var}\left\lbrack V \middle| z \right\rbrack}}} \\ {{{v(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}(z)}}}};} \end{matrix}$ a qα calculating step of calculating, via the computing device, a percentile qα of a portfolio value based on a first derivative 1′(z) of 1(z) calculated in the conditional expected value calculating step, a second derivative 1″(z) of 1(z) calculated in the conditional expected value calculating step, v(z) calculated in the conditional variance calculating step, and a first derivative v′(z) of v(z) calculated in the conditional variance calculating step using ${\Delta\;{q_{\alpha\; i}(z)}} = {V_{i}\frac{\partial\;}{\partial V_{i}}\left\{ {{- \frac{1}{2\;{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)} \right\}}$ ${\Delta\;{q_{\alpha}(z)}} = {\sum\limits_{i}{\Delta\; q_{\alpha\; i}}}$ ${q_{\alpha\; i} = {\left. \left( {{l_{i}(z)} + {\Delta\;{q_{\alpha\; i}(z)}}} \right) \middle| {}_{z = {N^{- 1}{({1 - \alpha})}}}q_{\alpha} \right. = {\sum\limits_{i}q_{\alpha\; i}}}};$ an E[V_(i)] calculating step of obtaining, via the computing device, an exposure value Vir of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, and obtaining a probability p_(ir) of transition from a rating i to the rating r stored in the storage device, and calculating an expected value E[V_(i)] of an exposure value of the obligor i based on the obtained V_(ir) and probability p_(ir) using ${{E\left\lbrack V_{i} \right\rbrack} = {\sum\limits_{r}{p_{ir}V_{ir}}}};\mspace{14mu}{and}$ a portfolio VaR calculating step of calculating, via the computing device, a portfolio VaR from the percentile q_(αi) and the expected value E[V_(i)] of the exposure value of the obligor I, wherein the rating threshold calculating unit includes: a ps_(sr) calculating step of calculating, via the computing device, a probability ps_(sr) of transition from the rating s to the rating r or lower based on the probability p_(sk), using ${p\; s_{sr}} = {\sum\limits_{k = 1}^{r}p_{sk}}$ a ps_(ir) calculating step of obtaining, via the computing device, a probability psi, of transition of the obligor i to the rating r or lower based on the probability ps_(sr) calculated in the ps_(sr) calculating step; and a θ_(ir) calculating step of calculating, via the computing device, the rating threshold θ_(ir) at which the rating of the obligor i becomes the rating r or lower based on the probability ps_(ir) obtained in the ps_(ir) calculating step using θ_(ir)=N⁻¹(ps_(ir)).
 7. A program product comprising a non-transitory computer usable medium having control logic stored thereon for causing a computer to function as: a z calculating unit calculating, based on a confidence level α stored in a storage device, a percentile of a risk factor z corresponding to a using z=N⁻¹(1−α); a rating threshold calculating unit obtaining a probability p_(sk) of transition from a rating s to a rating k stored in the storage device, and calculating a rating threshold θ_(ir) at which the rating of a obligor i becomes a rating r or lower based on the obtained p_(sk); a conditional probability calculating unit obtaining a constant a_(i) of each obligor i stored in the storage device, and calculating a probability ps_(ir)(z) that the rating of the obligor i becomes the rating r or lower under a condition that z is fixed based on the obtained a_(i), z calculated in the z calculating unit, and the θ_(ir), calculated in the rating threshold calculating unit using ${p\;{s_{ir}(z)}} = {N\left( \frac{\theta_{ir} - {\sqrt{a_{i}}z}}{\sqrt{1 - a_{i}}} \right)}$ where N is a cumulative probability function of a standard normal distribution; a w_(ir) calculating unit calculating a difference w_(ir) in exposure value between adjacent ratings of the obligor i based on an exposure value V_(ir) of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, using w_(ir)=V_(ir)−V_(i(r+i)); a conditional expected value calculating unit calculating a conditional expected value 1(z) of a portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating unit and ps_(ir)(z) calculated in the conditional probability calculating unit using $\begin{matrix} {{l_{i}(z)} = {\sum\limits_{r}{w_{ir}p\;{s_{ir}(z)}}}} \\ {{{l(z)} = {{E\left\lbrack V \middle| z \right\rbrack} = {\sum\limits_{i}{l_{i}(z)}}}};} \end{matrix}$ a conditional variance calculating unit calculating a conditional variance v(z) of the portfolio value V under a condition that z is fixed based on w_(ir) calculated in the w_(ir) calculating unit and ps_(ir)(z) calculated in the conditional probability calculating unit using $\begin{matrix} {{v_{i}(z)} = {V_{i}\frac{\partial\;}{\partial V_{i}}{{Var}\left\lbrack V \middle| z \right\rbrack}}} \\ {{{v(z)} = {\frac{1}{2}{\sum\limits_{i}{v_{i}(z)}}}};} \end{matrix}$ a qα calculating unit calculating a percentile qα of a portfolio value based on a first derivative 1′(z) of 1(z) calculated in the conditional expected value calculating unit, a second derivative 1″(z) of 1(z) calculated in the conditional expected value calculating unit, v(z) calculated in the conditional variance calculating unit, and a first derivative v′(z) of v(z) calculated in the conditional variance calculating unit using ${\Delta\;{q_{\alpha\; i}(z)}} = {V_{i}\frac{\partial\;}{\partial V_{i}}\left\{ {{- \frac{1}{2\;{l^{\prime}(z)}}}\left( {{v^{\prime}(z)} - {{v(z)}\left( {\frac{l^{''}(z)}{l^{\prime}(z)} + z} \right)}} \right)} \right\}}$ ${\Delta\;{q_{\alpha}(z)}} = {\sum\limits_{i}{\Delta\; q_{\alpha\; i}}}$ ${q_{\alpha\; i} = {\left. \left( {{l_{i}(z)} + {\Delta\;{q_{\alpha\; i}(z)}}} \right) \middle| {}_{z = {N^{- 1}{({1 - \alpha})}}}q_{\alpha} \right. = {\sum\limits_{i}q_{\alpha\; i}}}};$ an E[V_(i)] calculating unit obtaining an exposure value Vir of the obligor i when transition from a current rating to the rating r occurs, which is stored in the storage device, and obtaining a probability p_(ir) of transition from a rating i to the rating r stored in the storage device, and calculating an expected value E[V_(i)] of an exposure value of the obligor i based on the obtained V_(ir) and probability p_(ir) using ${{E\left\lbrack V_{i} \right\rbrack} = {\sum\limits_{r}{p_{ir}V_{ir}}}};\mspace{14mu}{and}$ a portfolio VaR calculating unit calculating a portfolio VaR from the percentile q_(αi) and the expected value E[V_(i)] of the exposure value of the obligor I, wherein the rating threshold calculating unit includes: a ps_(sr) calculating unit calculating a probability ps_(sr) of transition from the rating s to the rating r or lower based on the probability p_(sk), using ${p\; s_{sr}} = {\sum\limits_{k = 1}^{r}p_{sk}}$ a ps_(ir) calculating unit obtaining a probability psi, of transition of the obligor i to the rating r or lower based on the probability ps_(sr) calculated in the ps_(sr) calculating unit; and a θ_(ir) calculating unit calculating the rating threshold θ_(ir) at which the rating of the obligor i becomes the rating r or lower based on the probability ps_(ir) obtained in the ps_(ir) calculating unit using θ_(ir)=N⁻¹(ps_(ir)). 